平 谱 射 电 类 星体 B3 0307+380 的 15 GHz 射电 光 变 分 析 * 
HE, AB+ , eH 


( 云南 师范 大 学 物理 与 电子 信息 学 院 ， 云南 昆明 650500) 


摘要 : 欧文 斯 谷 射电 天 文 台 (OVRO) 的 长 时 间 监 测 数据 显示 平 谱 射 电 类 星体 B3 0307+380 
的 15 GHz 射电 光 变 曲线 中 表现 出 多 个 爆发 过 程 . 基于 双 指 数 函 数 ， 拟 合 了 12 个 爆发 过 程 ， 


得 到 了 每 个 爆发 过 程 的 光 变 幅度 和 上 升 、 下 降 时 标 ， 进 而 估算 了 亮 温 度 T, 和 光 变 多 普 勒 因 


Tó. 多 普 勒 因子 范围 为 (9.12 圭 0.58) ~ (35.38 土 3.34)， 平均 值 6, 为 18.65 土 7.11Csd). 


这 个 结果 表明 此 源 的 射电 辐射 存在 明显 的 队 束 效应 ， 文 持 主流 的 相对 论 性 喷 流 模型 ,利用 
Lomb-Scargle 周期 图 (LSP) 方 法 和 加 权 小 波 分 析 (WWZ) 方 法 ， 探讨 了 B3 0307+380 是 否 存 
在 周期 性 光 变 ， 发 现存 在 一 个 较 高 置信 和 度 〈( > 99% ) 的 约 244 天 的 周期 ， 并 简要 讨论 了 周 
期 可 能 的 存在 原因 . 
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中 图 分 类 号 : P162 ”文献 标识 码 : A ”文章 编号 : 


EZR AS (Blazar) 是 一 类 喷 流 与 视线 方向 近似 平行 的 活动 星系 核 ， 其 观测 特征 包括 宽 波段 
电磁 辐射 、 高 光度 、 高 偏振 、 视 超 光速 运动 等 ”， 耀 变 体 包括 平 谱 射 电 类 星体 (FESRQs) 与 蝎 
虎 型 天 体 (BL Lacertae Objects) 两 个 子 类 ， 两 者 的 差别 主要 体现 在 光谱 中 有 无 宽 发 射线 ". 
普遍 认为 ， 炮 变 体 在 儿 乎 整个 电磁 波段 都 表现 出 各 种 时 标 〔 从 几 分 钟 到 儿 年 量 级 不 等 ) 的 
光 变 ， 少 数 炮 变 体 表现 出 周期 性 光 变现 象 “”， 光 变 研 究 一 直 是 探索 炮 变 体内 部 辐射 过 程 和 


物理 机 制 的 最 有 效 手段 之 一 ， 由 于 次 变 体 噶 流 的 方向 与 观测 视线 之 间 的 夹 角 小 (0 <10")， 
炊 变 体 的 辐射 因 多 普 勒 效应 而 被 增强 ”"， 多 普 勒 因子 (5 ) 可 以 通过 喷 流 中 物质 流速 度 
CB) WAI CO) 这 两 个 本 征 参数 来 定义 ， 然 而 这 两 个 量 都 是 不 可 直接 观测 的 量 ， 因 
此 用 新 的 其 他 方法 来 估算 多 普 勒 因子 就 显得 很 必要 ， 在 各 种 估算 多 普 勒 因子 的 方法 中 ， 通 
过 射电 光 变 估算 多 普 勒 因子 (6, ) 的 方法 相对 比较 简单 ， 因 而 被 广泛 应 用 


B3 0307+380 (RA = 47.7078, Dec 238.2483 ) 是 一 个 被 探测 到 人 徊 马 射线 辐射 的 平 


谱 射 电 类 星体 "， 红 移 z = 0.816， 本文 利 用 欧文 斯 谷 射电 天 文 台 (Owens Valley Radio 
Observatory， 缩 写 为 OVRO) “的 15 GHz 长 时 间 监 测 数据 〈 观 测 时 长 约 11.4 年 ， 观 测 时 间 
从 2008 年 1 月 8 日 到 2019 年 5 月 25 日 ) ， 研 究 B30307+380 的 射电 光 变 ， 估 算 射 电 亮 温 
度 和 光 变 多 普 勒 因子 ， 并 考察 是 否 存 在 周期 性 光 变 .采用 的 宇宙 学 参数 为 


H, =73km:s™ :Mpc Q pu =0.27, Q =0.73. 


vacuum 


* 基 金 项 目 : 国家 自然 科学 基金 (11863007) 资 助 ; 云南 省 高 校 高 能 天 体 物 理 重点 实验 室 项 目 资助 . 
作者 简介 : 杨 星 ， 男 ， 硕 士 研究 生 ， 研 究 方向 : 活动 星系 核 . 
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1 数据 


利用 OVRO 的 官网 Chttp://www.astro.caltech.edu/ovroblazars/) ， 获 取 了 B3 03074380 
的 15 GHz 长 时 间 测 光 数 据 ， 去 除 少量 的 坏 点 后 ， 最 后 采用 的 测 光 数据 包括 664 个 点 ， 约 
MW Aig A MID 范围 约 为 54473 - 58628， 图 1 给 出 了 光 变 曲线 ， 横 坐标 为 MJD， 纵 坐标 为 以 


Jy 为 单位 的 15 GHz 流量 及 误差 (大 部 分 误差 很 小 ， 无 法 在 图 中 分 辨 ). 


2 多 普 勤 因子 估算 


基于 光 变 时 标 所 限定 的 尺度 ， 文 [13] 给 出 了 泡 变 体 亮 温度 的 计算 公式 : 


2 
T, = (4.5x10"°)AF 2D K (1) 
t (+z) 


sth, TAWE, AF 以 Jy HREM OGERE , t, DUROS I E, A DURS 


米 计 的 观测 波长 ，D 以 Mpc 计 的 光度 距离 . 
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图 1 B3 0307+380 的 15 GHz 光 变 曲线 
Fig.1 The 15 GHz light curve of B3 0307+380 


图 1 显示 B3 03074380 的 射电 光 变 曲线 中 表现 出 多 个 爆发 过 程 . 我 们 从 中 挑选 出 12 个 
包含 上 升 和 下 降 阶 段 的 完整 爆发 过 程 ， 采用 文 [14] 提 供 的 双 指 数 函数 公式 拟 合 这 12 个 爆发 
WHE. 具体 公式 如 下 : 


to—t t-t 


FQ)=F.+F(e" +e" J” 2) 


式 中 ， 丰 .表示 基底 流量 ，t 表示 峰 对 应 的 时 间 ， t, Alt, 分 别 表示 指数 上 升 时 标 和 下 降 时 


标 ， 五 ,衡量 爆发 的 幅度 . 


我 们 利用 基于 Python 语言 的 LMFIT 程序 包 进行 非 线 性 最 小 二 乘 拟 合 ， 调 用 的 函数 是 
leastsq0， 详 情 请 参考 https:W/lmfit.github.io/lmfittpy/， 拟 合 的 基本 原理 是 使 得 观测 流量 值 


N 
f t) 与 函数 推断 流量 值 严 (z,) 的 残 差 平方 和 SSR 达到 最 小 ，SSR = Y Lf) - FUL, 


i-l 


N 为 观测 点 个 数 . Levenberg-Marquardt (LM) 算 法 是 一 种 求解 非 线 性 模型 的 数值 算法 ， 常 


被 应 用 于 最 小 二 乘 曲线 拟 合 "5 . 我们 调用 的 函数 leastsq0 也 是 基于 这 种 算法 . 


12 个 爆发 过 程 的 拟 合 详情 和 参数 展示 在 图 2 和 表 1 中 .第 1 列 为 爆发 的 序号 ; 第 2 
列 为 爆发 对 应 的 约 化 儒 略 日 范围 ;第 3 列 为 观测 点 个 数 N ; 第 4 列 为 拟 合 所 得 的 约 化 最 


小 残 差 平方 和 启 ， | SSR 为 最 小 残 差 平方 和 ， 自 由 度 d.0.f = N-5; 第 5 


至 第 9 列 给 出 了 5 个 拟 合 参数 及 其 误差 ， 利 用 拟 合 参数 瓦 ，t Wr,» URAR (1) 计算 


了 每 个 上 升 或 下 降 阶段 对 应 的 亮 温 度 T,. 比 对 拟 合 公 式 (1) 与 文 [7] 中 的 公式 (2) ,我 


们 强调 在 计算 每 个 上 升 或 下 降 阶段 对 应 的 亮 温 度 T, 时 , AF 值 宜 取 拟 合 参数 五 的 一 半 . 


一 般 认 为 ， 粮 变 体 的 喷 流 成 分 对 应 的 亮 温度 克 不 能 超过 逆 康 普 顿 极限 温度 


T, =10"“K"““， 或 不 能 超过 平衡 亮 温 度 T = Sx10" K. Stokes T, 


约 分 布 在 (3.79 土 0.73)x10”~ (2.2140.63)x10°K ， 超 过 Ti 或 7 儿 个 量 级 ， 表 明 存 


著 的 多 普 勒 增 亮 效应 ， 类 似 文 [8]， 我 们 选取 也 , 作为 内 佘 亮 温度 .利用 如 下 公式 ， 
进一步 估算 了 对 应 的 光 变 多 普 勒 因子 : 


T, 
à, = e n (3) 


eq 


图 3 给 出 了 多 普 勒 因子 ó, 的 分 布 情况 .Oy 的 值 分 布 在 (9.12 土 0.58) ~ G5.3843.34) ， 
平均 值 为 18.65 土 7.11C(sd)， 中 值 约 为 18.14. 


表 1 爆发 过 程 的 拟 合 结果 
Table 1 The results of flare fitting 


MJD 
B f, F, F, t, f, 
No. range N 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
1 54696.4 - 
106 1.73E-03 54958.98+10.53 0.0140.05 1.5640.08 106.61+13.14 87.94+6.25 
55080.4 
2 55139.2 - 
555193 44 7.42E-04 §5228.7445.75 0.37+0.03 0.71+40.05 36.19+6.22 34.29+4.89 
3 55409.5 - 
41 3.91 E-03 55521.03+8.63 0.1140.09  1.82x0.17 47.2447.54 84.41412.07 
55665.8 
4 55666.8 - 
74 9.87E-04 55797.15x5.91 0.4140.01 0.73+40.03 60.11+5.28 40.344+4.84 
55969.9 
5 56060.9 - 
38 1.84E-03 56201.10+43.22 0.28+0.02 1.41+40.08 20.8622.36 64.36+6.19 
56325.9 
6 56330.3 - 
42 5.78E-04 56382.04+7.73 0.3840.01 — 0.8620.02 39.05+6.16 44.63+6.12 
56586.5 
7 56529.4 - 
61 6.44E-04 56704.68+3.50 0.3940.01  0.78+40.03 25.01+42.29 47.1043.75 
56859.5 
8 56941.3 - 
ease 2T 5.53E-04 57035.45+4.46 0.19x0.02 0.65+0.07 79.13412.44 19.22x3.28 
9 57200.6 - 
35 7.28E-04 57313.5743.28 0.1540.02 0.98+40.05 27.9642.71 132.94+9.36 
57531.7 
10 57597.8 - 
32 1.14E-03 57706.77+27.12 0.15+0.17 1.07+0.39 65.35+30.74 110.41+30.88 
57813.2 
11 58090.2 - 
22 5.61E-04 58157.86+8.50 0.20+0.03 1.22+0.06 52.55+7.67 41.70+6.46 
58309.5 
12 58276.7 - 
audies 6.59E-04 58449.49*7.22 0.2040.05 0.72+0.06 67.12+20.63 23.44+43.77 
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图 2 12 个 爆发 过 程 的 曲线 拟 合 
Fig.2 Curve fitting of twelve flares 
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图 3 光 变 多 普 勒 


因子 6 的 分 布 


Fig. 3 The distribution of variability Doppler factor Oy 


3 周期 分 析 


历史 文献 中 不 断 有 关于 炮 变 体 周期 性 光 变 的 报道 ”*”， 净 变 体 的 周 
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期 性 光 变 研究 ， 对 


理解 其 内 部 辐射 的 物理 机 制 和 辐射 区 的 几何 性 质 有 着 重要 的 意义 ， 目 前 ， 
常用 的 周期 分 析 方法 有 Lomb-Scargle 周期 图 (LSP) 方 法 和 加 权 小 波 分 析 方 法 (WWZ) ^ 
等 ， 我 们 也 利用 这 两 种 方法 探索 了 B3 0307+380 的 15 GHz 光 变 曲线 是 否 存在 周期 光 变 ， 结 


在 天 文学 上 比较 


果 列 在 图 4 F. 图 4 左 图 为 LSP 分 析 结 果 ， 图 中 我 们 按照 文 [23] 提 供 的 基于 红 噪 声 的 蒙特 卡 
洛 方法 给 出 了 9590. 99% 、99.7% 置 信 度 曲线 .结果 可 看 出 ， 存 在 一 个 约 244 士 11 天 的 
> 99% 短信 和 度 的 周期 ， 同 时 也 存在 一 个 约 465 土 40 天 的 >99% 置信 度 的 周期 , 后 者 大 约 为 


前 者 的 两 倍 ， 图 4 右 图 为 加 权 小 波 分 析 结 果 ， 


右 的 周期 ， 并 且 这 个 周期 基本 持续 ， 特 别 是 在 观测 时 间 的 中 段 . 


图 中 也 大 致 可 以 看 出 存在 一 个 大 约 250 天 左 


Normalized Power 


Time (days) 


图 4 LSP 周 
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期 分 析 结果 《〈 左 图 ) ; WWZ 周期 分 析 结 果 ( 右 图 ) 


Fig.4 The results of LSP (Cleft panel) and WWZ (right panel) period analysis 


总 结 与 讨论 


EN 


耀 变 体 光 变 曲 线 的 主要 特征 是 
相对 论 性 喷 流 的 聚 束 效应 是 导致 这 些 极端 观测 特征 出 现 的 原因 
(OVRO) 的 长 时 间 监 测 数据 显示 平 谱 射电 类 星体 B3 0307+380 的 15 GHz 射电 光 变 曲线 中 表 
我 们 利用 双 指 数 函 数 拟 合 了 12 个 爆发 过 程 ， 得 出 光 变 多 普 勒 因子 分 布 
在 (9.12 土 0.58) ~ B5.3843.34 范围 内 ， 平 均值 


现 出 多 个 爆发 过 程 


在 多 波段 呈现 出 大 幅 快速 光 变 且 存 如 


E 高 偏振 ， 通 常 认 为 
欧文 斯 谷 射电 天 文 台 


为 18.65 土 7.11C(sd). 这 表明 B3 


0307+380 的 射电 辐射 存在 显著 的 多 普 勒 增 亮 效应 ， 文 持 被 广泛 接受 的 相对 论 性 喷 流 模型 . 


粮 变 体 的 周期 光 变 可 能 与 钙 射 区 域 在 喷 流 
体 来 说 ， 由 于 螺旋 形 运动 导致 观测 视角 作 周 期 性 变化 ， 多 普 勒 
导致 多 普 勒 增 亮 的 辐射 流量 也 受 此 效应 的 调制 .图 5 给 出 了 多 普 勒 


间 的 变化 情况 ， 对 应 的 MID 值 取 拟 合 参数 加 及 其 误差 ， 多 普 勒 


P 治 着 螺旋 形 轨迹 作 相对 论 雪 
因子 也 发 生 相 应 的 变化 ， 从 而 
因子 在 12 个 爆发 过 程 之 


因子 Oy 取 上 升 和 下 降 阶 段 


运动 有 关 ” . A 


的 平均 值 及 其 误差 . 从 图 5 中 可 发 现 多 善 勒 因子 确 随时 间 呈 现 一 定 的 振荡 我们 猜测 B3 
0307+380 的 射电 周期 光 变 可 能 与 螺旋 喷 流 模型 有 关 .， 需要 指出 的 是 B3 0307+380 的 约 244 
天 的 射电 周期 置信 度 仅 为 大 于 2 水平， 还 有 待 进一步 确认 ， 另 外 ， 关 于 周 


喷 流 模型 仅 是 可 能 


原因， 还 有 很 多 其 他 的 机 制 值得 参考 


期 的 解释 ， 螺 旋 
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图 5 光 变 多 普 勒 因子 6, 的 变化 .每 个 点 对 应 的 MID 值 取 拟 合 参数 加 KERE, Oy 取 上 升 和 下 降 阶 段 对 


应 多 普 勒 因子 的 平均 值 及 其 误差 . 


Fig.5 The variability of Ó, . The MJD of each point adopts the fitting parameter fọ and its error, and Ó, 


adopts the average value of Doppler factors during the rise and decay processes and its error. 
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On the 15 GHz Radio Variability of  Flat-spectrum Radio Quasar B3 0307+380 
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Abstract: The long-term monitoring data of Owens Valley Radio Observatory (OVRO) show 
that FSRQ B3 03074380 has many flares in the 15 GHz light curve. Based on the double 


exponential form function, we have fit twelve flares, then derived the variability amplitude, the 


rise and decay time scales of each flare, and evaluated the brightness temperature (7,) and 
variability Doppler factor( ô, ). The range of ô, is from 9.12+0.58(sd) to 


35.38+3.34(sd), with the average of 18.65 - 7.11Csd?) . This result indicates that the radio 


emission is strongly beamed, and it is consistent with the well-known relativistic jet model. 
Making use of LSP and WWZ methods, we examined if there is a period in the radio light curve of 
B3 03074380. As a result, we find a period of about 244 days with a confidence level of 
>99% , and give the possible explanation about this period. 

Key words: Blazar; Brightness temperature; Doppler factor; Period 


